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Abstract 

A model is developed for the prediction of the seismo-acoustic noise spectrum 
in the microseism peak region (0.1 to 0.7 Hz). The model uses a theory devel- 
oped by Cato (J. Acoust. Soc. Am., 89 , 1096-1112 (1991)] for an infinite depth 
ocean in which the surface orbital motion caused by gravity waves may produce 
acoustic waves at twice the gravity wave frequency. Using directional wave spec- 
tra as inputs, acoustic source levels are computed and incorporated into a more 
realistic environment consisting of a horizontally stratified ocean with an elastic 
bottom. Noise predictions are made using directional wave spectra obtained from 
the SWADE surface buoys moored off the coast of Virginia and the SAFARI sound 
propagation code, with a bottom model derived using wave speeds measured in 
the EDGE deep seismic reflection survey. The predictions are analyzed for noise 
level variations with frequency, wave height, wind direction, and receiver depth. 
These predictions are compared to noise measurements made in ECONOMEX us- 
ing near-bottom receivers located close to the surface buoys. Good agreement is 
found between the predictions and observations under a variety of environmental 

conditions. 
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Chapter 1 



Introduction 



Researchers have long been aware of a peak in both the seismic ambient ground 
motion spectrum and the ocean bottom pressure spectrum which occurs with a 
period of about three to four seconds. This peak has historically been called the 
microseism peak. It has been noted that the frequency of this peak; seems to occur 
at twice the frequency of the peak in the surface gravity wave spectrum. Longuet- 
Higgins [lj was the first to develop a comprehensive theory which proposed the 
interaction of opposing surface waves, which closely approximate standing waves, 
as the source of the microseism peak. He showed that while the first-order pressure 
fluctuations caused by the gravity waves exhibit an exponential decay in depth, 
the second-order fluctuations caused by the nonlinear interaction of two opposing 
waves does not attenuate. Brekhovskikh [2] was the first to develop this idea into 
a prediction of noise in the ocean. He has been followed by others, most recently 
Kibblewhite and Wu [3], who all used a perturbation expansion as the solution 
to the wave equation with the resulting second order solution as the acoustic 
field. Cato [4] has developed a model which does not rely on a perturbation 
expansion, and does not require that a standing wave approximation be made. 
His theory predicts somewhat higher noise levels for a shallow receiver than those 
using the perturbation expansion, and a different directionality. In contrast to 
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these predictions, which have the wave interaction as the source mechanism, are 
those of Guo [5] who has proposed the wind turbulence acting on the surface 
directly as the source of noise in the microseism peak region. 

Previous measurements of noise in the deep ocean have been few, and the re- 
sults have been broadly consistent with all the above theories [6, 7, 8, 9]. What has 
been lacking is a direct accurate measurement of the source field believed to be 
causing the noise, either the wave field or the wind field. Without these quality 
source data, broad assumptions concerning the source field must be made, and 
the potential variance in the predictions based on these assumptions are typically 
greater than the differences in levels predicted by the different theories. Thus 
without the accurate measurement of the source field, it is difficult to evaluate 
the adequacy or inadequacy of particular theories. In the case of the surface 
wave field, the required measurements consist of directional surface wave spectra. 
This lack of high-quality, simultaneously measured, noise and surface wave data 
was the primary reason for the deployment of the ECONOMEX and SAMSON 
experiments [10,11]. An additional desirable feature of the ECONOMEX experi- 
ment was its long-term nature, allowing noise measurements under a wide variety 
of environmental conditions. The measurements consisted of near-seafloor water 
pressure and ground motion at various depths from ~100 meters to 2500 meters 
in the Atlantic Ocean. During the time the ECONOMEX instruments were de- 
ployed, instruments in the same area from the SWADE program were measuring 
surface parameters including directional wave spectra, wind speed, and other me- 
teorological quantities. These two coupled data sets provide researchers with the 
measurements needed to test the theories mentioned above. 

This work will examine the model developed by Cato [12] for an infinite depth 
ocean and will present the derivation of that model in detail. Using the SWADE 
directional wave spectra as inputs, we will estimate the resulting acoustic source 
level. We will then incorporate this source level into a more realistic environ- 
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ment consisting of a horizontally stratified ocean with an elastic bottom using the 
SAFARI program [13,14]. The geoacoustic model for our bottom comes in part 
from the EDGE deep seismic reflection study [15]. The resulting ambient noise 
predictions will then be compared to the measured noise from ECONOMEX in 
an effort to judge the adequacy of our model. 

In Chap. 2, we review the equations of motion and derive the inhomogeneous 
wave equation in terms of the Lighthill analogy [16]. The inhomogeneous wave 
equation is then solved in the presence of a sound speed and density discontinuity 
in Chap. 3. In Chap. 4, we make some simplifications to allow us to apply 
the theory to the motion of surface gravity waves. Additionally in this chapter 
we follow the work of Cato [12] in deriving an expression for the noise level. in 
terms of the directional wave spectrum, a problem for which we also introduce 
an elastic bottom. In Chap. 5 we describe the methods used to make predictions 
of noise levels using the equations in Chap. 4, the SWADE directional spectra, 
and the SAFARI program. Finally in Chap. 7 we present the ECONOMEX noise 
observations and compare them to our predictions which, for the most part, agree 
very well. Our conclusions and reco mm endations for future research are presented 
in Chap. 8. 
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Chapter 2 



Basic Equations 



In any consideration of sound generation, the most basic question to ask is “What 
is sound?” That is, how do we distinguish particle motions associated with acous- 
tic energy from that associated with other types of energy such as vortical energy 
(that due to turbulence) or thermal energy? In general, this is a very difficult 
question, for the equations governing the three types of motion are a set of cou- 
pled, inhomogeneous, non-linear partial differential equations, first derived in a 
linear form by Rayleigh [17], and later in a more general form by Doak [18]. Doak 
further showed that for small Stokes number the set decouples into three separate 
equations for the acoustical, vortical, and thermal type motions. For our model 
of sound generation, we will consider only this decoupled, small Stokes number 
approximation, and introduce any coupling between the three types of motion via 
external forcing, heat addition, or boundary conditions. The acoustic motion is 
then easily recognized as the small scale pressure perturbation associated only 
with the internal elastic properties of the medium. 

Let us consider the validity of the small Stokes number approximation in the 
case of low-frequency sound in the ocean. If we limit our interest to a maximum 
frequency of 50 Hz, then the maximum Stokes number (uu/c 2 ) is « 10 -7 , where w 
is the acoustic frequency, v is the medium viscosity, and c is the sound speed. We 
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are then justified in using the assumption of small Stokes number and can take 
advantage of the great mathematical and conceptual simplifications it allows. 

In our derivation of the sound generated by a moving fluid, we will follow the 
acoustic analogy method of Lighthill [16], but will use the pressure perturbation 
as our field variable in preference to fluctuations in mass density as suggested by 
Doak [18]. We will first derive an expression for the pressure perturbations in an 
ideal acoustic medium at rest, and then derive a similar expression for a real fluid 
which may contain variations in density, variations in sound speed, and which can 
have existing within it any type of motion. When we compare the two expressions, 
the difference will contain the effects of motion in the real fluid which generate 
sound. 

The exact equations of mass conservation and momentum can be written as 



dp_ d(pu,) 
dt dxi 



= 0 , 



( 2 . 1 ) 



d(pui) + djpujuj + a,j) _ = Q 



( 2 . 2 ) 



dt dxj 

where p is the mass density, u,- is the fluid velocity in the x t - direction, cr tJ is the 
stress tensor consisting of hydrostatic and viscous stresses, and /,• is the body 
force in the x,- direction. The indices i and j may take any value 1, 2, or 3, 
corresponding to the Cartesian coordinate axes, and a repeated index in a term 
indicates the term is to be summed over all values of the index. 

If we consider the “acoustic approximation”, p — po = c 0 2 (p — p 0 ) , where p is 
the fluid pressure, then Eq. (2.1) can be rewritten 

1 dp d(pui) 



2 p, H — ^ — 0. 

Co of dx{ 



(2.3) 



Doak has shown that this approximation is exact for an inviscid, non-heat con- 
ducting fluid, and valid for a Stokesian fluid to first order in the Stokes number. 
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Equation (2.2) can be applied to our ideal acoustic medium by linearizing it and 
neglecting viscous terms, leaving 

^ = 0. (2.4) 

dt dxi 

We can now combine Eqs. (2.3) and (2.4) by eliminating pu,- to yield 

d 2 p ,d 2 p 

aF~ c ^= °' (2 ' 5) 

This, then, is the familiar homogeneous wave equation governing the pressure 
perturbations in an ideal acoustic medium. 

Now consider the exact equation of momentum balance, Eq. (2.2). This can 
be expanded as 

d{pui) dp dipuiUj + Sij) 

ST + a7 i d7; + pf " (2 ' 6) 

where S is the viscous stress tensor, and the hydrostatic term has been explicitly 
stated. Next we again eliminate pu,- by using the exact equation of mass conser- 
vation, Eq. (2.1). By adding d 2 p/dx"l to both sides to allow direct comparison 
with Eq. (2.5), we arrive at 

d 2 P jd 2 p _ _ 2 d 2 (puiUj + Si,-) t d 2 (p-c 2 0 p) jd(pfi) , 0 ^ 

dF~ c «dxj c ° + ~~d? c “~- (2 ' 7) 

This is the inhomogeneous wave equation as derived by Doak [18]. It is iden- 
tical to that derived by Cato [4] with the exception of the final term, which he 
neglected. The pressure fluctuations in a real fluid are exactly those which would 
occur in a uniform acoustic fluid subject to the external stress system given by the 
right hand side of Eq. (2.7). As c% is assumed constant, not only the sound gener- 
ation but also the propagation effects of a real fluid are included in the equivalent 
stress system. It will be our task not only to solve Eq. (2.7) in an ocean geome- 
try with a realistic stress imposed, but also to separate out the sound generation 
terms from the propagation terms. 

It is useful at this point to consider the importance of the term c 2 d(pfi)/dxi, 
which was neglected by Cato in his theory of noise generated by surface orbital 
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motion [4]. In the ocean, this term reduces to c^g[dpjdz). We will consider 
motions with a periodic time dependence (p a e~' ut ) and a characteristic length 
scale L. Then we can form the ratio of the last to the first of the terms on the 
right hand side of Eq. (2.7) as ( gp/L)/(pu 2 /L 2 ). From Eq. (2.1), we can see 
\pU / L\ ~ \wp\ so we can rewrite our ratio as g/Lu> 2 . We can also see from Eq. (2.1) 
that the first two terms on the right hand side of Eq. (2.7) are of the same order. If 
we consider as our worst case a minimum frequency of 0.1 Hz, we can neglect the 
gravity dependent term if L 25m. If we were concerned with purely acoustic 
motion, the appropriate length scale would be the wavelength L = X = c/ f 
15,000 m. However the length scale appropriate to the equivalent sources resulting 
from surface gravity waves would be the wavelength obtained from the deep water 
low frequency surface gravity wave dispersion relation u 2 = 2%g/X which yields 
L = X cz 160 m. In both cases we can therefore neglect the last term in Eq. (2.7). 
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Chapter 3 



Solution to the Inhomogeneous 
Wave Equation 



The general solution to the inhomogeneous wave equation is well known (cf. Strat- 
ton [19]). For the source term in Eq. (2.7), it is given by 



PCM) ~Po = 

J_ r f d 2 (pu,u y + Sjj) d 2 (p/cg -p) ' dy 

47 tJ v [ dt/idyj dt 2 r r 



J_ r ri dp J_£r J_drdp 

"^47t J s [r dn r 2 dn ^ c^r dn < 

where r = |x-y|, h is the outward normal to the surface of integration, and 

the integrands axe evaluated at the retarded time r = t — r/c. The volume of 
integration above must include all possible regions where noise could be generated, 



dS( y), 



(3.1) 



and the surface integral will include the boundaries of these regions. We will choose 
our volume to be a cylinder of infinite height above the sea surface, and extending 



to the ocean bottom. The radius of this cylinder Y is assumed to be large enough 
to include any regions of motion which could contribute to p(x,t), and thus the 
sides of the cylinder do not contribute to the surface integral. We next follow 
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Cato (4] in modeling the sea surface density discontinuity as a Heaviside function 
so that the derivatives in Eq. (3.1) exist. The density is given by 

P = C) {Pw Pa) "h Pa 

where p w is the density in water, p a is the density in air, f is the value of x at the 
interface, and H(x — f) is the Heaviside function 

Ii(y) = 0, y > 0, 

H{y) = 1, y < 0. (3.2) 



We can then expand terms containing spatial derivatives of p as follows: 
d 2 (pUjUj) _ 

dyidyj 



q 2 

H— — (pty p a ) U t Uy + — — 



+2- 



dyidy, 
dH d 



dyidyj 



( Pw Pa)UjUj H“ (p«i Pa)UiUj H. 



dyidyj ^ r “'-'- J dy { dyj 
We can also combine Eqs. (2.1) and (2.2) to yield 

d 2 p _ d 2 (puiUj + <J it j) 



(3.3) 



(3.4) 



dt 2 dy { dyj 

which can be used with Eq. (3.3) to eliminate the derivatives of the Heaviside 
function in Eq. (3.1) resulting in 



4?r (p(x,t) - p 0 ) = 

f ^ 

J v r 



dy ^ H-z—z — (p w UiUj) + (1 - H)^-(paUiUj) + d2S<i 



dyidyj 



dyidyj dyidyj 



_f dy 
J v r 



I H - - ( p w UiUj ) + (1 - H) ^ (PaUiUj) + d - ^(4) 



dyidyj 



f T 1 dp 1 dr 1 dr dp 
J s [r dn r 2 dn ^ + CqT dn dt 



dyidyj 



dS( y). 



dy^yj dt 2 eg 



(3.5) 
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Now we can split the volume integral into one integral over the water volume 
and one integral over the air volume, noting that H = 1 in water and H = 0 in 
water, and after recombining terms to form d 2 p/dt 2 we find 



4tt(p(x, t) - po ) = 
f dy 

J v a r 



dy 

v w r 



+ 1 



d 2 . . d 2 . p . 

———(PvjUiUj + S{j) + (-j — p w ) 

ot/it/j m c 0 



d 2 p 

{p a UiUj + Si,) + — - Pa ) 



dyidyj 



f fl^p 1 dr dp 

/s[rdn r 2 dn^ CqT dn dt 



dt* K cl 



dS(y). 



(3.6) 



As before, the integrands must be evaluated at time r = t — r/c 

The next step is to apply the divergence theorem in a manner first described 
by Curie [20], and presented in detail by Cato [12]. This is not a straightforward 
process, as the terms to be evaluated depend on y both directly and also through 
the dependence on r. Curie stated and Cato demonstrated that 



r d 2 F i} - dy 


d 


r dFijdy 


( -d 


dFi, 1 


J vdyidyj r 


dx{J 


' v dy } - r 


1 vdy; 


. d yj r . 



dy 



f , dFjjdSjy) 

J s' d Vj r ’ 



(3.7) 



where /,• are the direction cosines of the outward normal of the surface enclosing 
V. We can then repeat this operation a second time to yield 



J v ay } - r oxjJ v r Js r 

We can now combine Eqs. (3.8) and (3.9) to achieve Curie’s result 
d 2 F ij dy _ 



(3.8) 



/ 



vdy { y y r 



oxiOXjJ v r ax{J s r J s oy } - r 



dy 



dF jj dS{y) 



(3.9) 
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We can also easily transform the surface integral in Eq. (3.1) as follows. First we 
resolve the normal derivatives into Cartesian coordinates which yields 



/ 



1 dp 1 dr 1 dr dp 
s [r dn r 2 dn P c%r dn dt 



dS( y) = 



/ li 

J s 



1 dp 1 dr 1 dr dp 
- + -—P + -J 



d5(y). 



■r dyi r 2 dy,- cgr dy,- dt J T 
By substituting dr /dyi = —dr/dxi and dp/dy,- = 5(p<5 tJ -)/y ; -, we have 

1 dr dp] 



f \l dp 1 dr 1 dr dp] _ 

J s'\r dy { r 2 dy P cgr dy,- dt J r y 



f d(p<? tJ ) d5(y) r | ! dr 1 dr dp 
J s * dy ; - r J s * [r da;,- cjj da:,- dt 



d5(y) 



Next we note that 






so that we can write 

d(p%) <fS(y) 



_ f i ZiE-iil zzMl — f i- z zL 
J s * dy ; - r J s * [r da:,- 



1 dr 1 dr dp 

p + - ^ 



dS(y) 



Cq dx{ dt 

d _ dS(y) 



f l d dSi < 

~ J s’ dx p8i ’ r 

Combining Eqs. (3.10), (3.11), and (3.12) yields 

f \l dp 1 dr 1 dr dp~\ . . 

/ |"o — I jT - ? H a - ^7 dS(y) = 

J s [r dn r 2 dn CQrdndt\ f 

J 5 ay ; * r J s ox{ r 



(3.10) 



( 3 - 11 ) 



(3.12) 



(3.13) 



Now we substitute Eqs. (3.9) and (3.13) into Eq. (3.5). After noting that the 
source motions at infinity do not contribute, we can separate the surface integrals 
enclosing the water volume into one at the air-sea interface designated 5,- and one 
at the ocean bottom designated St,, and can then combine the surface integral at 
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the bottom derived by applying the divergence theorem to V w with the surface 
integral of Eq. (3.13). This yields 



4tt ( p(x,t) - po ) = 



_d 2 

dx:dx 



*—[ (PvUiUj + Sij)— + / 

idxjJ Vy, r J Vy,ot 2 Cq r 



d 2 ,P 



dy 



, d ‘ f t , o \ d y , f fp \dy 

+ Wr>^ + ^ r + Jv.dt^4 Pa) r 



**/P 



,<*y 



- / s + s «) - [ s f i + s„) ™ 



<i5(y) 



dx, 

a 



s. l ’ {pMUi + + /,/' l4 ( '’“ UiU '' + y 

d f . . _ c ,dS(y) 

+ HtJ s> + S « + Mv)— — 

+ y" s ^> + S i> + Phi) 



dS(y) 



*S( y) 



(3.14) 



We can transform the y derivatives to time derivatives by noting from Eq. 
(2.2) that 

li 



— (pu.uy + Sq- + p£,y) 



- l i df (P u i)> 



and by defining Ap = (p„ — p a ), we can finally write 



47r(p(x,t) - P o) = 



_£ 2 

dx:dx 



idxj vS Pv,UiUi + Sii) r + Iv„dt 2 C 0 Pw) r 



d 2 ,p 



Ay 






3* ,P 









d5(y) 
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(3.15) 



We can now begin to interpret the meaning of Eq. (3.15) in terms of the sources 
of classical acoustics. The double spatial derivatives of the volume integrals are 
recognized as distributions of quadrupole sources in the volume V , and the single 
spatial derivatives of the volume integrals are recognized as dipole source distri- 
butions in the volume V. The spatial derivatives of the surface integrals at the 
interface are seen as dipole sources caused by the motion of the surface of the 
density discontinuity, and similarly the surface integrals of the time derivatives 
are seen as monopole source distributions caused by the motion of the surface of 
the density discontinuity. These interpretations are those given by Lighthill (16]. 

In order to interpret the terms involving the double time derivatives, we must 
make a further approximation to allow us to simplify Eq. (3.15), although it is 
one that is frequently made in classical acoustics. If we assume that the fluids axe 
isentropic, then the density and pressure are simply related by dp/dp = c 2 , where 
c is the local speed of sound [21]. Also, if entropy is conserved, the processes 
are reversible and the dissipative effects of viscosity can be ignored. (Note we 
must now treat the effects of attenuation due to absorption separately.) With no 
viscosity, 5 tJ = 0, or <7,y = p<5,y. We can thus write 

£ ( p_ _ \ £ Ap _ p. \ 

U ) W «v ’ 



or using Eq. (3.4) we have 





a 2 

dyidyj 



( PUiUj + Oij) . 



(3.16) 
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We can again apply the divergence theorem twice to Eq. (3.16) to yield 




(pUiUj + pS^ 

/ r . dy 

(puiUj + pbij) — 



pUi 



S(y) 



(3.17) 



From the above, we note that we can again identify both a quadrupole and a 
dipole volume source distribution, and both a dipole and a monopole surface 
distribution. Next we combine Eqs. (3.17) and (3.15) resulting in 



4?r(p(x, t) - p 0 ) = 



d 2 f \c 2 , . 1 rfy 



dy 

r 



+ a^Jv. + 

[(**4 - pA) u i u i + (4 - 4 ) p %] ^ 7 ^ 

, i [ , & \ ( 2 2\ i ds (y) 

+ clJs! i dt ^ PaCa)U ‘J r 



Co 

d 



f , cl, * . _ c ^ds{ y) 

■T" / h 2 (P* u i u i) + P 6 <: ) 

5 t eg r 

-/ 

J 5* Cq r 



(3.18) 



This result shows that the sound received at a point is the sum of apparent 
sources distributed throughout the source volume, at the air-sea interface, and at 
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the ocean bottom. The effects of the real ocean, such as refraction, are included 
in Eq. (3.18). The reflection from the ocean surface and bottom are also modeled 
in Eq. (3.18). Our next step will be to further simplify Eq. (3.18) to identify the 
real sources which transform other types of energy into acoustic energy, and to 
show how the sound generated by these real sources is modified by a homogeneous 
ocean. 
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Chapter 4 



Application to Orbital Motion in 
the Deep Ocean 



Equation (3.18) is quite general and could in theory be used to predict the acoustic 
pressure for a given source mechanism, but would require a knowledge of the 
particle motions u,- throughout the volume of the ocean and on the boundaries. 
In addition, the effects of a real fluid such as refraction are modeled as apparent 
sources, which further confuses the issue. We will now begin to make assumptions 
which will greatly simplify Eq. (3.18) to allow us to use it to predict infrasonic 
noise generated by the sea surface orbital motion. 

We first assume a plane wave sea surface elevation 

f = acos(/c -x* — at), (4*1) 



where f is the sea surface elevation, a is the amplitude, k is the wavenumber 
vector, x 1 = xii + x 2 j is the horizontal position vector, and a is the radian gravity 
wave frequency. The solution to the linear equations of motion yields a velocity 
potential <f> at depth x$ = z given by 



4 > = 



<racosh/c(z + d) 
k sinh Kd 



sin(fc • x' — ot), 



(4-2) 
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where d is the total water depth and k = |*c| [22]. Ignoring surface tension, this 
leads to the dispersion relationship 

a 2 = gK tanh Kd, (4-3) 

where g is the gravitational constant. If we restrict ourselves to water deep enough 
to satisfy Kd > 7 t / 2 , Eq. (4.3) reduces to the deep water surface gravity wave 
dispersion relation 

a 2 = Kg. (4.4) 

We then find the particle motions to be 

x\ = ae KZ ° sin(/c • xj, — at), i = 1,2, 

z = ae KZ ° cos(fc • Xq — at), (4.5) 

where the initial position of the particle is (x|), Zo). We can see that the particles 

describe a circular path, and hence the term orbital motion. Now we must de- 

termine the minimum water depth for which we may neglect the orbital motions 
at the bottom. We will also need to find the minimum water depth for which 
Eq. (4.4) holds. If we consider a minimum gravity wave frequency of 0.05 Hz, 
the deep water dispersion relationship is valid in water depths > 200m. To find 
the depth at which we can neglect orbital motion on the bottom, we consider the 
t~ KZ decay with depth and we find the orbital motions are 5% of the surface value 
at a water depth > 300m. We will therefore consider to be the effective source 
generation region to be <300 m. Outside this region we may neglect the volume 
integral of Eq. (3.18) and also neglect the orbital motions on the bottom. 

We can eliminate any apparent sources caused by refraction by considering 
the case of an isovelocity ocean, or c w = Co, and by considering the air also to be 
isovelocity, c a = q. We can also eliminate a number of terms in Eq. (3.18) by 
noting that c 2 w c 2 and p w c 2 w p a c\ . Thus, we neglect the air volume integral, 
and we assume the upper half-space to be a vacuum. We have now reduced the 
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motions on the bottom to those due only to the motion of the acoustic wave. 
Thus, in Eq. (3.18) we can eliminate from the bottom integral the p w UiUj and 
d(pu)/dt terms, since for an acoustic wave pu <C p. After all these simplifications, 
we are left with 



47r(p(x,f) -Po) = JL-f y 

~ j Ji + P%| 



One can now identify each term in Eq. (4.6) and relate it to a specific physical 
mechanism. The quadrupole volume integral involving pu,uy is identified as the 
volume contribution to the second order pressure fluctuation first proposed as 
the cause of microseisms by Longuet-Higgins (1] . The dipole surface integral 
involving pu,uy is identified as the sea surface contribution to this mechanism. 

r\ 

The monopole surface integral containing jfcpui is identified as the first order 
pressure effect which attenuates with depth. The apparent sources due to the 
surface and bottom reflections are contained in the two p<5 (J - terms . We can see 
in this case that the bottom does not influence the real source terms directly, and 
the real sources caused by the transformation of mechanical energy to acoustic 
energy can be evaluated first, and the bottom effects added later. 

Cato has shown [4] the quadrupole contribution to be < 2(t/ck times the dipole 
contribution in the far field, which is defined by r/c > 1. For a minimum wave 
frequency of 0.05 Hz, this is then valid at depths greater than 100 m below the 
region of effective volume generation, which we found to be no greater than 300 
m. Thus we are able to ignore the quadrupole term for receivers greater than 400 
m. Since we are mainly interested in the microseism peak region, we will also 
neglect the monopole term, although it will be significant at moderate depths for 
frequencies below the microseism peak. 



dS(y) 
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In finding the acoustic source level, we must make two major calculations. 
The first is to relate the source acoustic power density P(w,z) to the frequency- 
wavenumber spectrum of the source motions pu,uy. We must then relate this 
spectrum, say $(u;,k), to fi(oj which is the readily measured power spectrum of 
the sea surface elevation Both of these derivations have been performed by 
Cato [4,12], and we will present the important steps for the dipole sources below. 
We must then incorporate these source levels into an appropriate ocean model to 
find the noise power for an ocean receiver. 

We start by redefining some terms. Let 

pp(x,f) = J J,G(r)Wi,( y, r)dy, (4.7) 

where Wij = pu,uy, G(r ) = l/(47rr), and y is the source coordinate on the sea 
surface. Since we will perform the integration over the mean sea surface, we can 
immediately align the x s axis with the depth axis 2 and then li = I 2 = 0, l s = 1, 
and W{j = W{ = pu,-U 3 . We assume W{ to be temporally stationary and ergodic, 
and spatially homogeneous. Then the autocorrelation function of P£>(x,f) Ls 

Rp(t,t + c) = ((pd(x,0pd(x,« + ?)), (4.8) 

and since the process is ergodic 

Rp(t + c) = Pp(c) = lim -j- [ p D (x,t)p D (x,t + <;)dt, (4.9) 

T— *00 2,1 J -T 

and the power spectrum is 

r°° 1 r T 

Qd{u) = lim / — p D (x,t)p D (x,t + <;)dte' Wi d<; 

T— »oo J -00 Z1 J -T 

= r'iS, II w Hr [lk Is G{r)W, V- T)di , 
x \h- Is G W W '&' *wl ^ dc- (4-10) 
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Since y and y' are independent, we can define a separation parameter 1} in a cross 
correlation function of W, 



= Jim -L f T WMy.iJW', (?' + .), ( + <)*, (4.11) 

i — »oo J — T 

and we can further define the cross spectrum of 

A r°° 

$ii{v,u)= dq. (4.12) 

J—oo 

We also note that 

/ OO A]r 

(4.13) 

Finally we can combine Eqs. (4.10)-(4.13) and after making some substitu- 
tions, find that 

/ OO Air 

«„(w,k)fr l («,k,*)ir;(a;,k > *)— ( 4.14) 

00 (27TJ 

where 

Hi(oj,k,z) = (V ,k * f e -^r/c+k.t)^y \ ( 4 . 15 ) 

47T ox,- \ Js r J 

This completes the first step in determining the source levels from surface 
orbital motion. We have the source level Qd as an integral of the wavenumber- 
frequency spectrum of the orbital motions times the coupling factors H{Hf. Cato 
has calculated H analytically for an infinite ocean radius, and the values can be 
calculated numerically for a finite ocean radius [4]. Cato provides a thorough 
discussion of the behavior of the coupling factors [4] and has shown the spatial 
dependence of the coupling factors to be 

Hi oc exp t |k - x + z\J — A: 2 j j . (4-16) 

One can see that the spatial dependence is entirely oscillatory in the horizontal 
direction, but is oscillatory in the vertical only for k < u/c. For k > u>/c the 
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coupling factor magnitude has an exponential decay with depth. This will lead to 
an important approximation later in the development. 

Cato has also provided a convenient form for the coupling factors in terms of 
dimensionless quantities M = c o/kc, X = x?u>/c, X 0 = Rui/c and A — zw/c: 

z) = —icosaH 0 (u>,k,z)/M, 

H 2 {oJ,k,z) = — i since Ho(u,k,z)/M, 

Q 

H s (u,k,z) = —H 0 (oj,k,z), (4.17) 

where a is the angle between k and the t/i axis, and Hq is the monopole coupling 
factor 

1 f x o Xc _,VAr2+A2 

H a (w,k,z) = -f o -j= - = MX/M) dX, (4.18) 

Jo being the Bessel function of the first kind, order zero. 

We must now relate the frequency-wavenumber spectrum $,-j to the surface 
wave height power spectrum 17. We start by relating the spectrum of u,- to that 
of puiUj. Cato has shown [ 12 ] that 

$,/ lm (w,k)/p 2 = ^ ti (a;,k) * ^ Jm (w,k) + ^ fm (u;,k) * %(u;,k), (4.19) 



where * denotes convolution and (w, k) is the power spectrum of u,-. We can 
expand one of these terms in polar coordinates, with a being the gravity wave 
frequency, as 



[$«;jm(w,fc,a!)]i = P 2 [ f [ Vu{a,K,')) 

JkzzQ J 7=0 J o=—oo 



m — , # /cd/ccfadcr 

x — <T,K ,7 ) > 



(4.20) 



where k 1 = k — k, with 7 , 7 ' being the angles k, k 1 make with the yi axis. We can 
define & = w — < 7 , and can now identify o and & as the frequencies and k and k 1 
as the wavenumber vectors of the two interacting waves u t u ; -, and u> and k as the 
resulting acoustic frequency and wavenumber. We can relate the wavenumbers 
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and angles by 

kcosa = /ccos7 + k! cos 7', 

k sin a = Ksin') + k' sin')' . (4-21) 

Since a and k (and a 1 and «') axe associated with the orbital motion of surface 
gravity waves, they are uniquely related by Eq. (4.4). This is not true for oj and 
k. We can, however, use the behavior of the coupling factors previously described 
to find approximations for w and k in terms of a and k. 

We will again follow Cato [12] in noting the exponential decay of the coupling 
factors as k exceeds w/c, and we can choose a value (3 such that we will only 
consider a Fourier component to ensonify a region if k < fiio/c. From Cato’s 
examination of the coupling factors [4], we can choose (3 = 5 for a practical 
receiver depth. We can now write 

- Pw/c < |k| < fioj/c, (4.22) 

but |k| > |k| — |*:'| so we may use the deep water dispersion relationship to write 

— fioj/c < a 2 /g — (oj — o) 2 fg < flu/c, (4.23) 

or 

|u; — 2cr| < fig/c. (4-24) 

Thus, the acoustic frequency will be within Z% of twice the gravity wave frequency. 
We can say with reasonable accuracy that 

o» ~ 2cr ~ 2<7* . (4.25) 

We now have formally obtained the frequency doubling effect, where two gravity 
waves of the same frequency interact to generate an acoustic wave of twice that 
frequency. We can similarly show that k ~ k'. 

Next we define the two-dimensional wavenumber spectrum of u,- 

/ °° do 

( 4 - 26 ) 

-00 



27 



Since o and /c are related, we can define o as the specific value of a determined 
by Eq. (4.4). Also, since k and are not variables in the convolution integral Eq. 
(4.20), k' = |k — #c|, and we have the two relations of Eq. ( 4 . 21 ), we can find k 
and k' in terms of 7 . We will define the values of 0 and k thus determined as o 7 
and /c 7 . Now we will make the assumption that the two-dimensional wavenumber 
spectrum is separable into a one-dimensional spectrum times an angular spectrum, 



XilKl) = 2TTXii(/c )G(a,l)/K, 



(4.27) 



where the angular dependence satisfies 

[ G(c,rt)<h = 1 . 

JO 



(4.28) 



After substituting Eqs.(4.26) and (4.27) into Eq. (4.20), making the variable 
substitutions mentioned above, and also substituting d/c = (8k / do) do, we note 
that we can eliminate two of the three integrals in the convolution. Then 



($,yim(w,fc,a)]i = 2 tt p 2 f xu{K)X:m{K')G{o,^)G{o', 7 ') 

V-Y=0 



"1 

1 8 k . 



- d 7 . 



(4.29) 



'i.b') 3d, 

We arrive at a similar expression for the second term on the left hand side of Eq. 
(4.19) by substituting subscripts. We can write the sound pressure spectrum as 

r 2r 
'■ 7 = 



Qd M = [ 2irp 2 [ [x«i(/c)x/m(/c') 

J-00 J'l-O 



+ X.m(/c)Xl/(«')]G(o-,7)G(<7',7') 



X 



h* — 

do k' ij lm (2n) 2 



(4.30) 



Cato has shown that for significant ensonification of the noise field /c 7 , k and 
o can vary by only a small proportion of their values, and thus can be regarded 
as constants in the integration. Also, G(o, 7 ) varies relatively slowly, and k ~ /c', 
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so we can substitute G{o,^)G(o, 7 -f- tt) for G[o, 'i)G(o > , V). We can differentiate 
Eq. (4.5) to find 

X11 = cos 2 - 7 -X 33 , 

X22 = sin 2 7 * X 33 j 
X13 = cos 7 • X 33 , 

X23 = sin 7 * X33* (4.31) 

We further define the one-sided frequency spectrum of u 3 as fl(cr) and can write 

H(<r) = ^x»W, (4.32) 

and since u s = d$/dt, we have the frequency spectrum of the wave height in terms 
of the wavenumber spectrum 

n(<7) = (4 - 33) 

We can now see that 

Xu (*) = 9u ( 7 ) ^ (<?) !^o 2 <5 {/ , (4.34) 

where <7,/ can be found from Eq. (4.31) and <733 = 1. We note that 



9u ( 7 ) 533(7 + tt)^/ = ^.3 ( 7 ) ^ 3(7 + *) s u, 



(4.35) 



and so both terms of Eq. (4.30) are identical. We can then write the sound 
pressure in terms of fl(a), 

,a 4 dcr ft* f 2 * dadk 



_ , . ,a w .... , .... /■*« f 2r „ dadk „ , . 

Q D (w)=2p-—(l(a)I a Mj o f HiH.k—Sa, (4.36) 



where 



I c' it {°)=f (7) 1733 (7 + tt) ^(cr, 7) G(cr, 7 + 7r)d7. ( 4 . 37 ) 

Jo 

We can see that 7 aiI + I a22 = I a33 . 
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Next we separate out the a dependence of the Hi, perform the a integration 
of Eq. (4.36), and note that Hi = H 2 to find 

A 4 O A - 

Qd (w) = p 2 ^-£n 2 (a)I a3S (d) J Q [HiHi + 2 H z Hz) kdk , (4.38) 

where Hi = £T f (ty, z). Finally, we use the dispersion relationship to substitute 
k 1 = cr 12 / g, dcr/K = g/(2a ), and we use Eq. (4.25) to substitute 0 = 0 ' = w/2, 

and we have the result 
2 2 

Q d (u) = — wf2 2 (o;/2)J a3S (a;/2) f \kiH[ + 2H 3 H;)kdk. (4.39) 

4 Jo 

This gives us the received pressure for an infinite depth ocean in terms of known or 
measurable quantities. We can measure fi(cr) directly or use an empirical spectrum 
based on wind speed. We can similarly measure or imply from the wind speed 
G(cr, 7 ) and thus calculate I a33 (a). H{H ■ can be calculated for a given geometry, 
and p and g are physical constants. We thus have an equation which will allow us 
to predict the noise spectral levels for a receiver in an infinitely deep ocean. 

We must now turn to the problem of finding the noise levels in a more realistic 
ocean. We will consider the case of horizontally stratified media in the water 
column and in the bottom. Schmidt and Kuperman have shown that we can 
write the noise intensity generated by a horizontal distribution of homogeneous 
sources at depth z' as [23] 



P(u,z,z') = q 2 (u>,z')T2 l (u,z,z l ), 



(4.40) 



where q 2 (u>, z!) is the acoustic monopole source strength of the distributed sources 
of order m, and we can call T £ the bottom gain: 

T^(w,z,z') = R2m ^ J o \g{k,z,z')\[K 2 (z') - k 2 ]kdk. (4.41) 

In this equation K(z) = ui/c(z), and g is the depth-dependent Green’s function 
satisfying 

^ + [ K2 {z)-k 2 }g = -^-6(z-z'). (4.42) 
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We can eliminate the depth dependence in q(u>, z') by normalizing the source 
strength to that produced by the same source distribution in an infinitely deep 
ocean, <3(cu), where to first order q 2 ((j,z‘) = Q 2 (u)/16n(z') 2 . Now, provided z 1 is 
chosen small enough compared to the wavelength, Eq. (4.40) is approximately in- 
dependent of z 1 , and we have only to find the equivalent monopole source strength 
and evaluate T(u,z,z') to find P(cu,z). While the solution to Eq. (4.41) for an 
elastic bottom is not available analytically, there are numerical solutions available. 
We can easily relate Qd(w) to its equivalent monopole source strength by noting 
from Eq. (4.17) that 



HiHl + HtH; = HqHq/M 2 , 
oj d 

TT TT* TT TT* 

^ 3^3 — —-^—^ 0 ^ 0 - 



c dz 



(4.43) 



Thus the combined horizontal components are of vertical order m = 1 and the 
vertical component is of order m = 2. We can now write 



PdH = 



2 2 
P 9 



87T 



uCl 2 (u/2)I ass (cu/2) 



T 2 (u,z) H 0 HZ/M 2 kdk 
Jo 



+ T 2 (oj,z) [ k °-H 0 H*kdk 
Jo c 



(4.44) 



where the index on T indicates the order of the vertical model used in Eq. (4.41). 

We now have the full solution to the received pressure spectrum in a horizon- 
tally stratified ocean. We will next look at the specific techniques and data used 
to generate predictions of the noise pressure spectrum. 
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Chapter 5 



Prediction Techniques 



5.1 SWADE data 

The SWADE project (Surface Wave Dynamics Experiment) was an effort to char- 
acterize the sea surface using a variety of sensors and, at the same time, to measure 
other relevant environmental parameters (24). The project included several pitch 
and roll surface buoys, satellite radar backscatter measurements, SWATH ship ar- 
ray deployments, and aircraft overflights. The long-term deployment of the pitch 
and roll buoys is the element of the experiment of direct use to us in our effort to 
predict noise generated by surface orbital motion. 

SWADE and ECONOMEX instrument locations are shown in Fig. 1, and 
the overall experimental schematic is shown in Fig. 2. Deployment locations of 
SWADE instruments are listed in Table 5.1. 

The buoys provided one complete set of measurements including a spectral 



Instrument name 


Type 


Latitude 


Longitude 


Water depth 


DISCUS E 


Pitch and roll 


37° 20.0' TV 


73° 23.5W 


2670m 


DISCUS C 


Pitch and roll 


37° 32.1'TV 


74° 23.5W 


102m 


DISCUS N 


Pitch and roll 


38° 22.1'TV 


73° 38.9'W 


115m 



Table 5.1: Relevant SWADE instrument summary. 
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estimate each hour, allowing us to update our estimate of the noise each hour. 



The spectral data from the SWADE buoys was provided for frequencies from 0.03 
Hz to 0.34 Hz in 0.01 Hz bins, allowing acoustic predictions from 0.06 Hz to 
0.68 Hz. The SWADE buoys which were closest to the ECONOMEX instruments 
were the Discus E and Discus C buoys. For predicting the noise at the deep 
ECONOMEX site, the data of Discus E was used, while predictions of noise at 
the three shallow ECONOMEX instruments were based on the data of Discus C. 
The SWADE project had deployed a SPAR buoy near the shallow ECONOMEX 
site with better angular resolution than the Discus buoys, but it sunk prior to the 
ECONOMEX deployments. We should consider here the possible effect of using 
a surface buoy moored in 95 meters of water (Discus C) on the predictions of 
noise at the 450 meter and 790 meter OBS’s. We can note from the measured 
wave spectra that there is little energy below 0.1 Hz in the wave spectrum. If 
we calculate tanh(fcd) for the worst case of 0.1 Hz with d = 95 m we find our 
dispersion relationship is in error by less than 1%. Thus we would expect no 
correction need be made to our wave spectra measured at Discus C. 

Let us now consider the angular resolution available from pitch and roll mea- 
surements. Longuet-Higgins et al. [25] were the first to investigate the angular 
response of these buoys as follows. If the wavelengths of the surface motion are 
large with respect to the buoy diameter, the buoy tends to have the same motion 
and orientation as the surface. Then if we measure the vertical displacement and 
the two angles of pitch and roll, we will have three time series which represent 
the vertical displacement f, and its spatial derivatives d$/dxi and d$/dx 2 . Using 
the notation of the last chapter, we may represent the sea surface as a stochastic 
integral 



Then, since k = (/c cos 7, k sin 7), we can write our three time series denoted by 




(5.1) 
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6 , 6,6 as 



6 = 9?/* exp(t(K ' z — crt))dS 
Js 

6 = 9 lj i/c cos'] exp(t(K • x — at))dS 

6 = 9?/ t/c sin^y exp(t(#: • x — at))dS. 

Js 



( 5 . 2 ) 



Next we can form co-spectra Gij{a ) and quadrature spectra Qij{o) from the time 
series 6, 6, an< ^ we find 

r 2r 



Cn(a) = jf F(a, 7 )<f 7 , 

^(a) = f /c 2 cos 2 7^(0, 7) <£7, 

JO 

r 2 *- 

C'ssM = / k 2 sin 2 ')F(a, 7) <^7, 

JO 

r2x 

C 2S (a) = / /c 2 cos 7 sin 7-F(a, 7) <£7, 

r2r 

Qu{(r) = / /ccos7-F(a,7)d7, 

Jo 

f 2 * 

Qn(a) = / /csin7F(o,7)d7, 

JO 



( 5 . 3 ) 



where F(cr, 7) = n(o)G?(a,7) is the frequency-directional spectrum of the surface 
elevation. We can define the Fourier coefficients of F(a, 7) as 

1 /* 2 * 



where 



a n (<r) + t’6 n (a) = - f e n, ' , F(a, 7) <f7, 

7 T ^0 



n=oo 

= 2°o + S (“n cos(n-f) + 4 „sin(ti 7 )j. 



( 5 . 4 ) 



( 5 . 5 ) 



n= 1 



We can see that the right hand sides of Eq. ( 5 . 3 ) are related to these Fourier 
coefficients as follows: 



oo(o) = — C'u(a), 

7 T 

ai ( a ) = ^^ 12 ( a )’ 61 ^ 

a 2H = ^[^(a) - C'ssH), 62(a) = ^C 23 (a). ( 5 . 6 ) 
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The pitch and roll buoy then gives us the first five coefficients in the Fourier series 
describing the angular spectrum of the surface elevation at each of the frequencies 
for which we find the co-spectra and quadrature spectra of the time series. We 
must now use these five Fourier coefficients to find estimates of the frequency- 
directional spectra ft(er) and 

Our estimate of the wave power spectrum 

fl(<r) = f F((x,n) = nao{a) = C n (a), (5.7) 

Jo 

is obvious from inspection of Eq. (5.5), as well as our definition of the co-spectra. 

The best estimate of the directional spectrum is not as simple. An obvious 
choice would be to try the truncated sum 

G(a, 7) = — H — — («! cos 7 + 61 sin 7 + a 2 cos 27 + b 2 sin 27) . (5.8) 

2 7 T 7TOo 

This sum is actually a convolution of the true directional spectrum with a weight- 
ing function, and considerable smoothing results in the estimated spectrum. Cal- 
culations of I a3S made using this type of directional spectrum estimate from simu- 
lated directional spectra are typically in error by a factors of 10 s . Other weighted 
averages of the first five Fourier coefficients can be made, but they too produce a 
much smoothed estimate. 

Several investigators have fit empirical curves to measured directional spectra. 
Longuet-Higgins has suggested the wave directional spectra fit the form [25] 

G(o, 7) a | cos 25 (7/2) | (5.9) 

where the spreading parameter s is a function of frequency and wind speed. Kib- 
blewhite and Wu used an empirical relationship to find s based on the wind speed 
and frequency, and then calculated I azz analytically [3]. One could also match 
the measured first five Fourier coefficients to the first five Fourier coefficients of 
the empirical spectrum cos 2 ' 7/2 to estimate the parameter s, as suggested by 
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Longuet-Higgins [25], and hence calculate I a33 . However, the data used to de- 
velop the empirical formula are generally taken under conditions of steady wind 
speed and direction. In our field data, the wind speed and direction can vary 
significantly, giving rise to wave fields with different directionality, and in gen- 
eral a broader directional spectrum than that predicted by the empirical formula. 
Therefore, we would expect predictions of I a33 calculated from spectra derived 
from empirical formulas to be lower than the true value under variable meteoro- 
logical conditions. In particular, if the true directional spectrum is bimodal (two 
peaks corresponding to two wave fields generated by winds in different directions) 
there can be significant energy in opposing wave directions G(^)G(^ + 7 r) which is 
not predicted by the cosine power curve. Donelan ef al. have suggested a better 
fit to the data is found in a sech 2 /? / 7 distribution [26] with (5 being the spreading 
parameter, but estimates of I a33 based on estimates of (5 are also too small. 

Another approach would be to use a data adaptive spectral estimation tech- 
nique such as the maximum likelihood method (MLM) or the maximum entropy 
method (MEM). The method used here to estimate the directional spectrum given 
the first five Fourier coefficients is the MEM. This method produces a spectal es- 
timate which retains the first five Fourier coefficients and estimates the remaining 
coefficients based on the first five. We will follow Lygre and Krogstad [27] in 
developing an algorithm to make this estimate. 

We will define a function with a Fourier series (suppressing the dependence on 
<j) on the interval (— 7r,7r) as 

1 00 

Dfr) = TZ £ C. ‘o = 1, c_„ = <. (5.10) 

n= — oo 

The entropy of D is defined by 

H(D) = 1 r log (£>(7))*Y, (5.11) 

J — T 

and it has been shown by Burg [28] that the function maximizing H(D) subject 
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to the constraint that the coefficients c n equal some known c * for k < N is 

1 



D(l) = 



0] 



2tt |1 - (f> le-*^ <f> N e~ iN '*D\ 2 ’ 

where <f > i • • • ^ and 0* are obtained from the Yule- Walker equations 



(5.12) 



and 
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(5.13) 



*2 = 1 - M 






w 



(5.14) 



In our case we have N = 2 with ci = (a! + ibi) / a 0 and c 2 = (<*2 + *^ 2 ) /«o- We can 
now solve this system of equations to find <f> n in terms of c n : 



<f>i = (ci - c 2 cj)/(l - |c x | 2 ), 

<t> 2 = C 2 — Ci^i, (5.15) 



and finally we can substitute these into Eq. (5.12) to find our directional spectral 
estimate 



Dh) = 



1 — ^iCi — <^2 C 2 

2 tt |1 — — <^ 2 e 2t ^| 2 * 



(5.16) 



We now have an estimate of the angular distribution of the wave energy at each 
frequency which reproduces the first five Fourier coefficients which produced it, 
and uses these and the Yule- Walker equations to extrapolate the remaining co- 
efficients. This technique has been shown by Lygre and Krogstad [27] to give a 
much more peaked distribution over the MLM technique, and to resolve a bimodal 
wavefield. 

Figure 3 shows the estimated spectra from the various methods for a simulated 
bimodal spectral input, and Fig. 4 shows the estimates of directional spectra for 
the 0.16 Hz bin from Discus E on January 27 at 1200. This datum was chosen as 



37 



illistrative because it occurs after a shift in wind direction, and one could expect 
the true spectrum to be bimodal. One can see for the estimates from the simulated 
data that the MEM estimate provides a better representation of the structure of 
the actual spectrum. One can also note in the comparison of actual data a bimodal 
structure is evident in the MEM estimate only. 



5.2 Coupling factors and bottom gain 

While the coupling factors / HiH-kdk can be evaluated analytically for an ocean 
of infinite radius, [4] they predict infinite noise when the horizontal components 
are considered. We can solve this by adding a relaxation mechanism and thereby 
attenuation in the water column, or we can use a finite effective radius with some 
physical basis such as the ocean basin radius, storm radius (if applicable), or such. 
The coupling factors must then be evaluated either approximately analytically or 
numerically. The attenuation due to absorption at the frequencies of interest is 
believed to be too small to limit the noise on ocean basin scales, and so we will 
assume an effective radius. Calculations show that at distant ranges the horizontal 
coupling factors vary roughly a s y/R, and thus doubling the effective radius has 
only a minor effect on the overall source strength Qd{v). 

Based on the distance to shore of the ECONOMEX deployments and on trans- 
mission loss studies in this frequency range, an effective radius of 100 kilometers 
was chosen. The magnitude of the coupling factors was then calculated numeri- 
cally for receiver depths of 450 meters and 2500 meters. These results axe shown 
in Fig. 5. 

The final step in making our predictions is the calculation of the bottom gain, 
T £, for the horizontal and vertical components. The most computationally inten- 
sive part of this task is the calculation of the depth dependent Green’s functions, 
g(k, z, z'), for which we use the fast field approach. The fast field approach can be 
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generalized to a fully elastic elastic media, but requires that the wave equation be 
separable in depth, that is to say range-independent. We will use the fast field ap- 
proach, that realizing errors that may ensue due to range dependent bathymetry 
and media. 

The tool used to solve the full wave problem in this work is the SAFARI set of 
programs (Seismo-Acoustic Fast field Algorithm for Range Independent environ- 
ments) developed by Schmidt and Jensen [13,14]. The basic solution technique in 
SAFARI is to represent the field in each of a series of homogeneous layers by the 
Hankel transforms of the unknown potentials satisfying the homogeneous wave 
equation. The boundary conditions at each layer interface yield a set of local 
equations involving the unknown potentials of the adjacent layers. These local 
equations are collected into a global matrix which can then be solved to yield all 
the unknown potentials simultaneously. Solutions are determined efficiently by 
implementing modem numerical techniques. An additional advantage relating to 
this work is the inclusion of the integration of Eq. (4.41) as an option in the code. 

The environmental model used in SAFARI requires a number of parameters. 
The user must specify for each layer the compressional wave speed, shear wave 
speed, compressional and shear attenuations, and density. The uppermost and 
lowermost layers are taken to be semi-infinite half-spaces. Due to the long wave- 
lengths involved at the very low frequencies of our predictions, the environmental 
model should be as accurate as possible fairly deep into the bottom. It is fortu- 
nate that compressional wave speed data from a deep seismic reflection study of 
the U.S. mid-Atlantic continental margin was made available prior to publication. 
The EDGE seismic experiment [15] involved recording seismic profiles to 16 sec- 
onds off the Virginia coast in the same region as the ECONOMEX and SWADE 
experiments. Using the compressional wave speed data, the figures and equations 
of Hamilton [29] and discussions with other investigators working in the area of 
geoacoustic modeling [30,31], geoacoustic models were developed for the bottom 
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layer depth (m) C p (m/s) 


C, (m/s) 


Ip (dB jX) 


1. (dB/ A) 


p (g/ cm 3 ) 


vacuum 


* 


0 


0 


0 


0 


0 


fluid 


0 


1500 


0 


.015 


0 


1 


elastic 


2500 


1756 


200 


.5 


1 


1.8 


elastic 


3165 


2193 


540 


.3 


.5 


2.0 


elastic 


3678 


2495 


895 


.1 


.3 


2.2 


elastic 


5011 


2648 


1314 


.1 


.3 


2.4 


elastic 


5276 


3817 


1900 


.1 


.3 


2.7 


elastic 


7231 


4042 


2011 


.1 


.3 


2.8 


elastic 


8896 


6380 


3180 


.05 


.1 


3.0 



Table 5.2: Deep site environmental model. C Pi , and 7 Pi , are compressional and 
shear wave speeds and attenuations, respectively; p is the density. 



layer 


depth (m) 


C p (m/s) 


C, (m/s) 


% ( dB/X ) 


7, (dB/X) 


p (g/ cm 3 ) 


vacuum 


* 


0 


0 


0 


0 


0 


fluid 


0 


1500 


0 


.015 


0 


1 


elastic 


450 


1799 


200 


.5 


1 


1.8 


elastic 


970 


2174 


540 


.3 


.5 


2.0 


elastic 


1396 


2683 


895 


.1 


.3 


2.2 


elastic 


2441 


3441 


1720 


.1 


.3 


2.4 


elastic 


3332 


4398 


1739 


.1 


.3 


2.5 


elastic 


4938 


5872 


2933 


.1 


.3 


2.7 


elastic 


6582 


6159 


3059 


.05 


.1 


2.9 


elastic 


8227 


6380 


3170 


.05 


.1 


3.0 



Table 5.3: Shallow site environmental model. C Pt „ and 7 Pi< are compressional and 
shear wave speeds and attenuations, respectively; p is the density. 



at locations corresponding to water depths of 450 meters and 2500 meters. The 
models thus developed are given in tables 5.2 and 5.3. Depth profiles of com- 
pressional and shear wave speeds at the two sites are given in Fig. 6. SAFARI 
calculations of for both m = 1 and m = 2 for a receiver at 450 meters and 
2500 meters are given in Fig. 7. 

We can obtain a qualitative understanding of the relative importance of the 
different propagation mechanisms by looking at Fig. 8, contour plots of the magni- 
tude of the wavenumber integrands plotted against the inverse of the phase speed, 
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or “slowness” k/u, and frequency. One can see for the deep site the normal modes 
play the predominant role in propagating the surface noise to the deep receiver, 
as only for the lowest frequencies do we get a contribution from waves with phase 
speeds less than 1500 m/s. In the shallow case there is a significant contribution 
at all the frequencies of interest from the lower phase speeds, indicating the im- 
portance of interface waves. These waves are propagating horizontally but suffer 
an exponential decay in the vertical, and thus they are excited only in the shallow 
case. 

We can combine the effects of the coupling factors and the bottom gain, along 
with the constants in Eq. (4.44) to find the temporally invariant part of the 
solution, and we can write our prediction as the one-sided sound pressure spectrum 
level (dB re lfiPat/Hz) 



SL{f) = 20 log H(//2) + 10 log /«„(/) + 101og£(/), (5.17) 

where f is the acoustic frequency / = w/(2x) and B(f) is this time invariant part: 
10log2?(/) = 197 + 10 log J/ 

H$H$k dk )] • (5-18) 

The quantity 101ogi?(/) is shown in Fig. 9 for both the deep and shallow sites. 



+ T Uf)f 

Jo 



Tl{f) [ k ° H Q H*/M 2 kdk 
J o 
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Chapter 6 



Predictions for Receivers at 450 
meters and 2500 meters 



6.1 Variation of spectral level with frequency 
(spectral shape) 

Since the wave height power spectrum appears in Eq. (4.19) as fl 2 (/ /2), we might 
expect the shape of the acoustic spectrum to be related to the shape of the wave 
height spectrum at double the wave frequency, but we will see this shape is mod- 
ified by a number of factors. The wave height spectrum generally shows a very 
steep rise to a spectral peak followed by a somewhat gentler (oc o’ -4 ) slope at 
frequencies above the peak [26], This peak frequency is generally characterized as 
being inversely proportional to wind speed, and thus is usually lower in frequency 
at higher wave heights. The peak tends to be quite narrow, normally occurring 
at frequencies of about 0.1 to 0.2 Hz, although in a newly developing wave field 
it can be higher. This tends to give rise to an acoustic spectrum with a peak in 
the 0.2 to 0.4 Hz range. 

This overall shape will of course be modified by the effects of the directional 
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spectrum, the coupling factors, and the bottom gain as functions of frequency. 
Effects of the coupling factors are easily seen in Fig. 5 and tend to emphasize the 
higher frequencies, which will tend to mitigate the slope of the peak in the acoustic 
spectrum at frequencies above the peak. The spectral slope of the coupling factors 
is oc u> 2 . The bottom contribution to the spectral shape is more complicated, as 
seen in Fig. 7. In the shallow case, a peak in the bottom gain appears at about 
0.18 Hz, which sharply drops to a low at about 0.28 Hz, and follows with a rise 
above 0.28 Hz with a slope of about w 2 . This will tend to flatten the acoustic 
spectrum if the peak in the wave spectrum occurs above 0.1 Hz (as it almost 
always does). The peak in the bottom gain at 0.18 Hz is normally well overcome 
by the sharp drop in the wave height spectrum at wave frequencies below 0.1 Hz. 
In the deep case, the peak in the bottom gain occurs at about 0.24 Hz, and is more 
likely to have an effect on the shape of the acoustic spectrum. In most cases, the 
frequency of this peak in the bottom gain falls slightly below twice the frequency 
in the wave height spectrum, thus serving to broaden the peak in the acoustic 
spectrum. Above 0.3 Hz the bottom gain is relatively flat in the deep case. Thus 
in the deep location the bottom will tend to enhance the peak in the acoustic 
spectrum, and at frequencies above the peak the down slope should be greater at 
the deep location than in the shallower case. When we look at the overall transfer 
function B(f) shown in Fig. 9, we see for the shallow site a post peak slope of 
oc u s and for the shallow case a slope of oc tu 5 . 

By far the greatest effect on the shape of the noise spectrum next to the wave 
height power spectrum is the behavior of J aj3 . If one examines the behavior of 
the empirical models, it is clear that the value of the spreading parameters s or 
/? vary such as to reduce the value of I a3S at the frequency of the spectral peak. 
Mitsuyasu et al. [32] have proposed a model for frequencies above the spectral 
peak in which s varies as s oc ( U /ga)~ 2 5 , where U is the wind speed. In Donelan’s 
model, /? is dependent on cr/cr p only, with a p being the frequency at the peak [26], 
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Variation of I a33 versus frequency for this model is shown in Fig- 10. One can see 
that for either of these empirical models the directional spectrum is narrower near 
the peak frequencies, and broader away from the peaks, and thus the spreading 
integral I a33 will be a minimum near the peak. The effects of the wave height power 
spectrum and the spreading integral will oppose one another, and again tend to 
flatten the spectrum. We can see in Figs. 11- 26 the estimates of 10logP(/), 
n(tr/2n), 20log(fl(//2)), and 101og/ ass (/) for several illustrative examples for a 
2500 m receiver, and in Figs. 27 - 30 the same quantities for a 450 m receiver. One 
can see from Figs. 11- 30 that the wave height power input can vary by as much as 
60 dB in a given spectrum over the frequency range, whereas the spreading input 
varies over a much narrower range of up to 20 dB in a given spectrum. This will 
cause the wave height power spectrum to dominate, giving rise to a spectral peak, 
although one much diminished from the peak in the generating wave spectrum. 

In s umm ary we find the shape of the acoustic spectrum will in general resemble 
that of the generating wave spectrum in that their will be a sharp rise to a spectral 
peak, but in the acoustic spectrum the peak will be of a lower magnitude and 
broader than that present in the wave spectrum. 



6.2 Variation of spectral level with wave height 

The variation in the predicted noise with wave height will be the result of the 
combined effects of the wave height spectrum squared and the spreading integral. 
We saw in the last section that the two effects opposed one another, with the 
variation in the wave height spectrum dominating, giving rise to a spectral peak. 
When we look at predicted spectra from different times corresponding to different 
meteorological conditions, we find the variation in the wave height spectra is again 
greater than that of the spreading integral. We would thus predict an increase in 
noise level with an increase in overall wave height or sea state, with the increase 
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in wave height again partially offset by the decrease in the spreading. 

To allow us to see the effects of increasing wave height on our predictions, 
it is best to study a period of relatively constant wind direction with increasing 
wind speed and corresponding wave height. This will minimize the effects of 
other conditions which could effect the wave directional spectrum, such as wind 
direction versus fetch direction, sharp changes in wind direction, etc. We can see 
from Fig. 31 that we have such a time available to us at Discus E from 0300 to 
1200 on January 21, 1991. During that time, wind speed rises from about 5 m/s 
to about 12 m/s with a rise in significant wave height of 0.25 m/s to 1.9 m/s. 
Some predicted noise spectra for this period are included as Figs. 11-14. One can 
see in the noise prediction at 0300 a peak forming at 0.6 Hz with a level of 121 dB. 
In subsequent predictions, this peak moves lower in frequency, to 0.45 Hz at 0600, 
0.38 Hz at 0900 and finally 0.30 Hz at 1200. The peak level also increases from 
121 dB to 125 dB, 130 dB, and finally 133 dB at 1200. It is interesting to note 
the wave height power at the peak for this period varies from -13 dB to 12 dB, a 
range of 25 dB, while the value of the spreading integral at the peak varies from 
-18 dB to -25 dB, for a range of -7 dB. This clearly shows the mitigating effect 
of the spreading integral in limiting the noise power at the peak under moderate 
conditions. 

It is also interesting to consider the predictions in the case of severe weather 
conditions. The highest significant wave height recorded for which ECONOMEX 
data is available was 5.9 meters at Discus E in the late afternoon of March 4. 
Again the wind direction was relatively steady. The meteorological data for this 
period is shown in Fig. 32, and the noise predictions for 0400 to 1600 on that 
day are seen as Figs. 15-21. The peak noise level prediction for this time is 155 
dB, with a wave height input of 35 dB and a spreading integral value of -21 dB. 
Investigation of our MEM estimates of I a33 shows that it reaches a minimum of 
about -25 dB at the peak frequency under conditions of moderate wave height (~ 3 



45 



meters), with no further decrease with increasing wave height above that level. 
This means our estimated directional distribution is not becoming increasingly 
narrow at wave heights above ~ 3m. This is in disagreement with the observed 
behavior of the directional spectra of wave height at high sea state [26]. We 
might, therefore, expect our predictions to overestimate the noise at high sea 
states. The increase in predicted noise level at wave heights above this threshold 
will be directly proportional to the increase in the wave height power squared. 

At frequencies above the peak, the variation in wave height power with wave 
height is much less than that at the peak. One can see that the variation of the 
average wave height power on January 27 from 0300 to 1200 in the 0.5 Hz to 
0.68 Hz was only from ~ -20 dB to ~ -15 dB. This is typical of most of the data 
analyzed in this frequency range. Even under the extreme conditions of March 4 
the wave height power in this frequency range does not increase much above this 
level. Similarly the spreading input varies little in this range of frequencies for a 
steady wind direction, varying from ~ -5 dB to ~ 0 dB. We would therefore predict 
very little variation in the noise level in this frequency range for an increasing sea 
state. 



6.3 Variation of spectral level with changing wind 
direction 

One of the more interesting studies we can make involves examining the varia- 
tion of noise level under conditions of constant wind speed but changing wind 
direction. A theory which proposes that wind turbulence is the direct cause of 
the acoustic noise would predict very little variation in noise levels under these 
conditions. Thus the predictions we make here may be useful in judging the ac- 
tual contribution from the different mechanisms. Under our theory, we intuitively 
expect the noise to increase from the increased spread in the directional spectrum 
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as the new wave field is developed, but we also realize the original wave field is 
diminishing, due to the loss of wind forcing and the combined effects of dissipation 
and non-linear interactions with the newly developing wave field. 

We can see these combined effects in the data from Discus E in the early hours 
of February 23, 1991. As we can see in Fig. 33 the wind direction veers sharply 
by about 60 degrees at 0200 while the wind speed stays relatively constant at 8 
to 11 m/s. The response in the predicted noise spectra is shown in Figs. 22-26. 
One can see the noise prediction at the peak is fairly constant at about 135 dB 
from 0200 through 0400, then begins to rise until it reaches 140 dB at 0600. The 
wave height power input during this period actually drops from 15 dB to 8 dB 
at 0400, and it rises back to about 15 dB at 0600. The spreading input at the 
peak increases from about -23 dB at 0200 to -17 dB at 0400, and then remains 
at about this value through 0600. Thus the increase in predicted noise is due to 
the increased spread in the directional spectrum, not an increase in wave height. 
Another effect of the wind shift on the predicted acoustic spectrum is that the 
peak tends to broaden. This is seen most clearly at 0400 in Fig. 24, when the 
peak in the developing wave field and the peak in the pre-existing wave field are 
both of similar magnitudes but different frequencies. This behavior is typical of 
the wind shifts analyzed. 



6.4 Variation of spectral level with receiver depth 

The variation in predicted spectral levels with receiver depth can be attributed 
to three factors. First, for a receiver in very shallow water, the gravity wave 
dispersion relationship will begin to depart from Eq. (4.4) and the bottom will 
begin to have an effect on the wave height spectrum. For our receivers at 2500 m 
and 450 m, this is not the case. The two remaining factors, the coupling factor 
differences and the bottom gain differences, will cause variation in our expected 
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noise levels at these two depths. 

Referring to Fig. 5, we can see that the difference in the coupling factor gain 
between the two locations varies from about 15 dB at 0.06 hz to about 7 dB at 0.68 
Hz, with the shallow site having the higher value. The bottom gain is again more 
complicated, with the shallow bottom gain being higher at frequencies less than 
0.2 Hz and greater than 0.5 Hz, and the deep bottom gain being higher between 
about 0.2 Hz and 0.4 Hz. The combined effects are visible in Fig. 9, where we can 
see the differences in the energetic part of the spectrum are quite minimal. Given 
the same wave height directional spectrum input, we expect the deep case to yield 
higher noise levels of about 6 dB at 0.25 Hz, and we expect the shallow case to 
yield higher noise levels of about 5 dB at 0.68 Hz. Due to the limited amount, of 
data analyzed and the variance of the estimated spectra, it is difficult to see this 
variation in the predicted spectra between the two depths. We will therefore not 
try to judge the success of the predictions in this area. 
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Chapter 7 



Experimental Noise 
Measurements 



7.1 Experimental description 

ECONOMEX (Environmentally Controlled Oceanfloor Noise Monitoring Experi- 
ment) [10] was designed to provide a long-term, high quality seismo-acoustic noise 
data set which could be coupled to the surface wave and meteorological data of 
the SWADE experiment. The instrumentation consisted of six Office of Naval 
Research (ONR) ocean bottom seismometers (OBS’s) and two, one vertical and 
one horizontal, 75 meter six element hydrophone arrays. The instruments were off 
the Virginia coast in January 1991, recovered in February 1991 for maintenance 
, and redeployed from February through early April 1991. Precise instrument 
locations and deployment dates are listed in table 7.1. 

The ONR OBS instruments deployed consisted of a three-component geophone 
for measuring ground motion in the 0.07 to 80 Hz range, a Cox- Webb differential 
pressure gauge (DPG) for measuring long period pressure signals in the water 
column, and in the original deployment, an OAS hydrophone [34]. In the later 
deployment, these hydrophones were removed to improve instrument reliability. 
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Instrument 


Lat. (N) 


Lon. (W) 


Depth 


Deployment dates 


Vert, array 


37° 24.7' 


73° 26.8' 


2573m 


Jan 25-Feb 6; Feb 22-Apr 5 


Hor. array 


37° 24.7' 


73° 26.8' 


2573m 


Jan 25-Feb 6; Feb 22-Apr 5 


OBS 56 


37° 24.7' 


73° 26.8' 


2548m 


Jan 25-Feb 6; Feb 22-Apr 5 


OBS 58 


37° 26.4' 


73° 31.4' 


2417m 


Jan 24-Feb 7; Feb 22-Apr 5 


OBS 61 


37° 23.8' 


73° 24.4' 


2600m 


Jan 25-Feb 6; Feb 21-Apr 4 


OBS 62 


37° 33.2' 


74° 14.1' 


769m 


Jan 11-Feb 6; Feb 22-Apr 5 


OBS 63 


37° 34.1' 


74° 16.5' 


443m 


Jan 10-Feb 6; Feb 22-Apr 5 


OBS 51 


37° 35 .$ 


|74° 21.3' j 


95m 


Jan 10-Feb 6; Feb 22-Apr 5 



Table 7.1: ECONOMEX instrument summary. Positions and depths listed are for 
the second leg of the experiment; those for the first leg differ only slightly. Also 
note OBS frames 61 and 62 exchanged positions between the first and second legs, 
although the instruments on them were exchanged also such that the instruments 
remained deployed in the same locations. 



The sensors were connected via preamplifiers to an acquisition package consisting 
of a pre-whitening and anti-aliasing filter, a gain-ranging amplifier to improve 
dynamic range, and an analog to digital converter. The combined filter response 
is shown in Fig. 34. The acquisition package fed a recording package consisting of a 
RAM buffer and an optical disc recording system capable of storing 400 megabytes 
of data. The typical OBS’s were programmed for continuous 8 Hz recording, with 
the anti-aliasing filter set to 2 Hz; however OBS 56 was set to record at 128 Hz 
with its anti-aliasing filter set to 40 Hz. 

The 75 meter horizontal and vertical arrays each consisted of six OAS hy- 
drophones at 15 meter separation. The hydrophone signals were preamplified by 
a low-noise, wide-range preamplifier and sent to acquisition and recording pack- 
ages identical to those of the OBS’s, with 128 Hz sampling and the anti-aliasing 
filter set to 40 Hz. In the second deployment the bottom three hydrophones of 
the vertical array were not included due to a cable malfunction. The array cable 
jacket included loose ended fiber strands to reduce strumming noise. 
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7.2 Data selection and processing 

Following instrument recovery, the ECONOMEX data were transcribed from the 
optical disc to magnetic tape. While the data has not yet been transcribed from 
the binary machine format to a standard format for further dissemination, it 
is possible to read the binary format and produce ASCII files for limited time 
periods. The present work has concentrated on using the differential pressure 
gauge (DPG) data at the 8 Hz sampling rate to minimize the data processing 
involved while still adequately sampling the frequency band of interest. The DPG 
data were used instead of the geophone data because of current uncertainties in 
the geophone response. Since the model developed is only valid for receivers at 
depths greater than 400 meters, it was decided to analyze selected data from one 
deep DPG (OBS 58 at 2417 m) and the 450 m DPG (OBS 63). This allowed the 
maximum depth variation comparison given the instrument deployment depths. 

Estimated spectra were generated by removing the mean and any linear trend 
from 64 second segments of the time series, and then averaging 512 point fast 
Fourier transforms of 34 minute sections of data using a Hanning window. The 
response of the pre-whitening and anti-aliasing filter was then removed. The 
bandwidth of 64 seconds results in a frequency resolution of 0.0156 Hz, which 
is consistent with the spectral resolution of our predictions, which is 0.02 Hz. 
An example of a full spectral range observation in shown in Fig. 35. Since the 
predictions are limited to a 0.06 to 0.68 Hz band, further observed spectra shown 
in this work are limited to the same band for clarity. For the most part, data were 
analyzed which corresponded to the same time periods for which predictions were 
made in Chap. 6. 
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7.3 Observed results and comparison with pre- 



dictions 

7.3.1 Spectral shape and overall spectral noise levels 

In general, there is good overall agreement between the predicted spectral levels 
at most frequencies in the band of interest. Fig. 36 shows an example of one such 
case at the 2500 meter site, and Fig. 37 shows similar results at the 450 meter site. 
While the levels at a given frequency may differ between predicted and observed 
by up to 5 dB, there is an overall correspondence between the two. It is interesting 
to note in many of the comparisons between observed and predicted spectra that 
small peaks exist in the observed levels which are present in the predictions at 
the same frequenies but with differnt magnitudes. These two figures represent 
examples of the best agreement between predicted and observed spectra. 

More typical of the level of agreement are Figs. 38-43, where we can see close 
agreement at frequencies around the spectral peak, but differences away from the 
peak of up to 7 dB in the deep case, and up to 17 dB in the shallow case. Here 
the overall shape is correctly predicted, but the peak is broader or narrower in 
the observed spectra, giving rise to large differences in the high slope region. In 
the predictions from the shallow site, there is a tendency to predict levels that are 
too high in the band 0.1 Hz to 0.2 Hz. Since this band corresponds to the peak in 
the bottom response in the shallow case, errors in the bottom model may account 
for this difference. 

7.3.2 Variation with wave height 

The observations show an increase in noise level with wave height, with good 
agreement between observations and predictions at moderate (~l-3 meters) wave 
heights, as seen in Figs. 39-43. At very low predicted noise levels, which corre- 
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spond to times of low wave height, there is a disparity between predictions and 
observations. An example of this is the developing wave field of the early hours 
of January 27 at Discus E. As one can see in Figs. 44- 46, the agreement is good 
at frequencies corresponding to the peak in the predictions as the wave field de- 
velops, but the low levels predicted away from the peak are not confirmed by 
observation. By 1200 the wave field is developed, and Fig. 39 shows the good 
agreement between the prediction and observation at this time. A possible cause 
of the error in the predictions at low wave heights is the assumption of spatial ho- 
mogeneity of the source wave field. Under very low local wave height conditions, 
it is possible for a much stronger wave field at some distance to dominate the noise 
field, thus making the predictions made from the local wave field very much in 
error. Another possible cause of the differences under low wave height conditions 
is the existence of another source mechanism generating acoustic energy, whose 
noise is normally dominated by that caused by orbital motion. Under low source 
strength conditions for the orbital motion noise, this assumed source may now 
dominate, giving rise to the errors noted above. The predictions in general are in 
reasonably close agreement when measured significant wave height is above ~ 1 
meter. For the period of the ECONOMEX data, roughly 77 percent of the wave 
height measurements are above this threshold. 

We can see from the data of March 4, as seen in Fig. 47-50, the predictions again 
begin to deviate at very high wave height conditions which correspond to strong 
winds at relatively constant direction, with the predictions being higher than the 
observations. The disagreement at frequencies corresponding to the peak is up to 
10 dB under these conditions. There is, however, the same general trend in the 
observed noise data as exists in the predictions, that of higher levels at higher 
wave heights and wind speeds. This would tend to indicate the wave directional 
spectral estimate provided by the MEM technique is overestimating the spread at 
these high wave heights. The empirical models, which are based on strong, steady 
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winds, predict less spreading under these conditions than do the MEM estimates. 
Thus, under these conditions, the empirical model estimate for the spreading may 
be more accurate. The observations match the predictions fairly well in the region 
of frequencies above the peak under high wave height conditions. 

7.3.3 Variation with changes in wind direction 

The predictions agree quite well with the observed data under conditions of chang- 
ing wind direction. We can again consider the wind shift of February 23 at 0200, 
depicted in Fig. 33, as typical. The observed and predicted noise levels for sub- 
sequent times are shown in Figs. 51- 55. The observed noise spectrum changes 
with time roughly as the predicted spectrum does. Of particular note is the broad 
peak in the observed spectrum at 0400, and the higher levels seen at the peak 
from 0500-0700. 

It is also interesting to note the predictions and observations during a time 
of highly variable conditions such as those of the afternoon of January 27, again 
depicted in Fig. 31. The observations and predictions are seen in Figs. 36,56,57. 
Again we see relatively close agreement between predictions and observations, 
with the increased spreading adding to the noise as the wave height diminishes to 
keep the overall noise level fairly constant. 



54 



Chapter 8 



Conclusions and 
Recommendations for Future 
Research 



8.1 Conclusions 

The most important conclusion one can draw from this work is the apparent im- 
portance of the orbital motion contribution to the total noise level in the frequency 
band 0.1 to 0.7 Hz. The close correspondence between the predicted and observed 
spectra under a wide range of conditions is strong evidence that the true noise 
generation mechanisms in this band axe dominated by the interaction of opposing 
surface gravity waves. The fact that the predictions hold under the conditions of 
changing wind direction but constant wind speed would tend to negate the impor- 
tance of the direct input of the wind turbulence on the sea surface as an important 
sound generation mechanism under most conditions. There is a possibility that 
some other source mechanism or mechanisms contribute significantly to the noise 
spectrum in this band under conditions of low wave height. 

The second conclusion one may draw from this work is the importance of the 
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measured directional spectrum in determining the overall noise level. Any model 
based on empirical relationships using wind speed as their input would have to also 
include the growth and decay of the wave field based on changing wind direction. 
The models which predict directional wave spectra from a knowledge of wind speed 
alone would be doomed to failure in the general case due to their assumption of 
a fully developed wave field and the changing nature of the true wave field. 

Another important point brought out by this work is the importance of the 
propagation mechanisms in determining the overall levels and shape of the acous- 
tic noise spectrum in this frequency band, which is consistant with the work 
of Schmidt and Kuperman [23]. It is imperative that one take into account the 
propagation if one wishes to compare source levels between two different locations. 
Otherwise differences in the bottom contribution could cloud important correla- 
tions or lead to incorrect conclusions. A key element in research of this type is 
the geoacoustic bottom model, and the availability of measured wave speeds deep 
into the bottom in the present work was extremely fortunate. 

8.2 Recommendations for future research 

As the predictions of acoustic noise depend so critically on the angular spread of 
the wave energy, a fruitful line of investigation would be to obtain directional wave 
height spectra of greater angular resolution in future experiments. It is unfortu- 
nate that the Spar buoy, the sensor with the greatest directional capability in the 
SWADE project, was lost prior to the ECONOMEX deployments. The increased 
directional resolution could have been used not only in estimating spectra at the 
location of the buoy, but also could have helped in determining the error in the 
estimate of the directional spectra made when we use the MEM technique at other 
locations. In lieu of higher resolution directional wave spectra, an attempt could 
be made to use some of the other SWADE data, such as directional spectral es- 
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timates made by shipborne arrays or radar backscatter, to improve the estimates 
from the pitch and roll buoys. 

Another investigation which should provide interesting results would be the 
correlation of large amounts of the ECONOMEX data with the various SWADE 
parameters such as wind speed and wave height, and perhaps to automate the gen- 
eration of acoustic predictions from the SWADE directional spectra. The present 
research only scratched the surface of the data available from this long term ex- 
periment. It would be useful to find the amount of long term agreement, and even 
more interesting to find other periods of disagreement between the predictions 
made with the model developed herein and ECONOMEX observations. 

The final recommendation for research in this area would be the investigation 
of the range dependent aspects of the problem. There is range dependence in 
the both the source mechanism and the propagation mechanisms, which perhaps 
gives rise to some of the disagreement between our predictions and observations. 
The shallow site, in particular, would benefit from a consideration of its range 
dependent bathymetry. Buckingham [33] has found the noise in a wedged-shaped 
ocean with pressure-release boundaries to closely approximate the noise field in 
the range independent case, but further work is needed to extend his work to a 
more general bottom. As well, the investigation of the spatial variation of the 
wave height spectrum should prove feasible once the entire SWADE project data 
is collected and correlated. This could then be used to more accurately predict 
the noise at the ECONOMEX sites. 
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Latitude 



ECONOMEX and SWADE instrument locations 




Longitude 



Figure 1: Approximate locations of the ECONOMEX and relevant SWADE instru- 
ments. denotes pitch and roll buoy location, “x” denotes location of 
one ONR OBS and two hydrophone arrays, and “o” denotes the location 
of one ONR OBS. 
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Figure 2: Schematic arrangement of ECONOMEX and SWADE deployments. 
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Figure 3: Estimates of the directional spectrum of (7(7) = A r (cos 8 (7/2)-bcos 10 (7/2+ 
7r/4)) , where 7 is the azimuthal angle, and N is chosen to normalize the 
spectrum. The MEM estimate was made using the Lygre-Krogstad al- 
gorithm mentioned in the text. The MLM estimate was made using an 
algorithm by Lacoss [35]. The LCS estimate refers to the weighted aver- 
age of the first five Fourier coefficients suggested by Longuet-Higgins et 
al. [25] The empirical estimate was made by estimating the parameter 
s from the first five Fourier coefficients and using the empirical formula 
<7(7) = W cos 2 * (7/2). 
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Magnitude of G 



Directional spectral estimates for Discus E at 0.16 Hz, 1/27/91 1200 




Figure 4: Estimated spectra derived from the 0.16 Hz data bin from the Discus E 
buoy on January 27, 1991 at 1200 hours. The estimates were made using 
the methods discussed in the caption to Fig. 3. 
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Coupling factor integrals; 2500 m and 450 m receivers 




Figure 5: Plots of lOlog + H 2 H^)kdk and lOlog/* 0 H 3 H^kdk for re- 
ceiver depths of 450 meters and 2500 meters, ko was taken as 5w/c, and 
an effective radius of 100 kilometers was used. 
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Figure 6: Profiles of compressional and shear wave speeds, Cp and Cs, used in the en- 
vironmental model. Compressional wave speeds are based on unpublished 
data from the EDGE deep seismic reflection survey, and the shear wave 
speeds were derived from the compressional wave speeds, as described in 
the text. Figure 6a is the deep site and Fig. 6b is the shallow site. 
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Bottom Gain 




Figure 7: Plots of T % from data generated by SAFARI for the assumed bottom 
model at 450 m and 2500 m receiver depths for m = 1 and m— 2. 
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Figure 8: Contours of the integrands of Eq. (4.41), the depth-dependent Green’s 
functions g{k,z } z'), versus horizontal slowness 1 /c — k/o; and frequency 
for a) the deep site (2500 m) and b) the shallow site (450 m). 
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Figure 9: Plots of 10 log B(f) for the deep and shallow sites. 
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Spreading integral versus frequency using Donelan’s empirical model 




Figure 10: Plot of 10 log J a33 versus normalized frequency // /o using Donelan’s em- 
pirical model for the directional wave spectrum [26] where / 0 is double 
the frequency at the peak in the wave height spectrum. 
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20 log Wave Ht. Power (dB) Pred. Noise Power (dB) 






Figure 11: Predicted noise level and inputs for a 2500 meter receiver on 1/27/91 at 
0300. a) Noise spectral level 101ogP(/) in dB re /iPa 2 /Hz. b) Wave 
height power spectrum fi(a/ 2?r) in m 2 /Hz. c) Wave height acoustic input 
201ogfi(//2) in dB re m 2 /Hz. d) Spreading integral level 101og/ a33 (/) 
in dB re Hz -1 . 
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Figure 12: Predicted noise level and inputs for a 2500 m receiver on 1/27/91 at 0600. 

a) Noise spectral level 10 log P(f) in dB re /xPa 2 /Hz. b) Wave height power 
spectrum Cl(a/2ir) in m 2 /Hz. c) Wave height acoustic input 20logn(//2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log / ass (/) in dB re Hz -1 . 
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Figure 13: Predicted noise level and inputs for a 2500 m receiver on 1/27/91 at 0900. 

a) Noise spectral level 10 log P[f) in dB re /zPa 2 /Hz. b) Wave height power 
spectrum n(cr/2rr) in m 2 /Hz. c) Wave height acoustic input 201ogn(//2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log /„„(/) in dB re Hz -1 . 
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Figure 14: Predicted noise level and inputs for a 2500 m receiver on 1/27/91 at 1200. 

a) Noise spectral level 10 log P(f) in dB re /iPa 2 /Hz. b) Wave height power 
spectrum n(<r/27r) in m 2 /Hz. c) Wave height acoustic input 20 log Q(// 2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log /<*„(/) in dB re Hz -1 . 
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Figure 15: Predicted noise level and inputs for a 2500 m receiver on 3/04/91 at 0400. 

a) Noise spectral level 10 log P(f) in dB re /xPa 2 /Hz. b) Wave height power 
spectrum Q(<7/2tt) in m 2 /Hz. c) Wave height acoustic input 20logQ(//2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log / aj3 (/) in dB re Hz” 1 . 
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Figure 16: Predicted noise level and inputs for a 2500 m receiver on 3/04/91 at 0600. 

a) Noise spectral level 10 log P(f) in dB re /xPa 2 /Hz. b) Wave height power 
spectrum Cl(a/2n) in m 2 /Hz. c) Wave height acoustic input 20 log ft(//2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log in dB re Hz -1 . 
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Figure 17: Predicted noise level and inputs for a 2500 m receiver on 3/04/91 at 0800. 

a) Noise spectral level 10 log P(f) in dB re /xPa 2 /Hz. b) Wave height power 
spectrum Q(a/2n) in m 2 /Hz. c) Wave height acoustic input 201ogfi(//2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log /«„(/) in dB re Hz” 1 . 
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Figure 18: Predicted noise level and inputs for a 2500 m receiver on 3/04/91 at 1000. 

a) Noise spectral level 10 log P(f) in dB re/iPa 2 /Hz. b) Wave height power 
spectrum fi(cr/27r) in m 2 /Hz. c) Wave height acoustic input 20log0(//2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log /*,,(/) in dB re Hz“ x . 
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Figure 19: Predicted noise level and inputs for a 2500 m receiver on 3/04/91 at 1200. 

a) Noise spectral level 10 log P(f) in dB re /iPa 2 /Hz. b) Wave height power 
spectrum ft(cr/27r) in m 2 /Hz. c) Wave height acoustic input 20logft(//2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log (/) in dB re Hz -1 . 
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Figure 20: Predicted noise level and inputs for a 2500 m receiver on 3/04/91 at 1400. 

a) Noise spectral level 10 log P( f) in dB re /xPa 2 /Hz. b) Wave height power 
spectrum Q(cr/27r) in m 2 /Hz. c) Wave height acoustic input 20 log Q(//2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log /<*,,(/) in dB re Hz -1 . 
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Figure 21: Predicted noise level and inputs for a 2500 m receiver on 3/04/91 at 1600. 

a) Noise spectral level 10 log P(f) in dB re /xPa 2 /Hz. b) Wave height power 
spectrum ft(<r/27r) in m 2 /Hz. c) Wave height acoustic input 20logfi(//2) 
in dB re m 2 /Hz. d) Spreading integral level 10log/ aj3 (/) in dB re Hz" 1 . 
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Figure 22: Predicted noise level and inputs for a 2500 m receiver on 2/23/91 at 0200. 

a) Noise spectral level 10 log P( f) in dB re /iPa 2 /Hz. b) Wave height power 
spectrum Q((r/2ir) in m 2 /Hz. c) Wave height acoustic input 20 log ft(//2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log (/) in dB re Hz -1 . 
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Figure 23: Predicted noise level and inputs for a 2500 m receiver on 2/23/91 at 0300. 

a) Noise spectral level 10 log P( f) in dB re /zPa 2 /Hz. b) Wave height power 
spectrum n(cr/27r) in m 2 /Hz. c) Wave height acoustic input 20logH(//2) 
in dB re m 2 /Hz. d) Spreading integral level lOlog J as3 (/) in dB re Hz” 1 . 
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Figure 24: Predicted noise level and inputs for a 2500 m receiver on 2/23/91 at 0400. 

a) Noise spectral level 10 log P(f) in dB re fi Pa 2 /Hz. b) Wave height power 
spectrum Q(<r/27r) in m 2 /Hz. c) Wave height acoustic input 20 log fi(//2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log I azz (/) in dB re Hz -1 . 
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Figure 25: Predicted noise level and inputs for a 2500 m receiver on 2/23/91 at 0500. 

a) Noise spectral level 10 log P(f) in dB re jxPa 2 /Hz. b) Wave height power 
spectrum Cl{a/2n) in m 2 /Hz. c) Wave height acoustic input 20 log fi(//2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log /„„(/) in dB re Hz -1 . 
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Figure 26: Predicted noise level and inputs for a 2500 m receiver on 2/23/91 at 0600. 

a) Noise spectral level 10 log P[f) in dB re /i Pa 2 /Hz. b) Wave height power 
spectrum Q(cr/2n) in m 2 /Hz. c) Wave height acoustic input 20 log H(// 2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log I aii (/) in dB re Hz” 1 . 
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Figure 27: Predicted noise level and inputs for a 450 m receiver on 2/23/91 at 1000. a) 
Noise spectral level 101ogP(/) in dB re /iPa 2 /Hz. b) Wave height power 
spectrum Q(a/2ir) in m 2 /Hz. c) Wave height acoustic input 201ogfi(//2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log /<*„(/) in dB re Hz -1 . 
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Figure 28: Predicted noise level and inputs for a 450 m receiver on 2/23/91 at 1400. a) 
Noise spectral level 101ogP(/) in dB re fi Pa 2 /Hz. b) Wave height power 
spectrum 0(cr/27r) in m 2 /Hz. c) Wave height acoustic input 20 log fl(//2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log / a33 (/) in dB re Hz” 1 . 
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Figure 29: Predicted noise level and inputs for a 450 m receiver on 2/27/91 at 1500. a) 
Noise spectral level 10logP(/) in dB re /zPa 2 /Hz. b) Wave height power 
spectrum fi(cr/27r) in m 2 /Hz. c) Wave height acoustic input 201ogQ(//2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log /<*„(/) in dB re Hz” 1 . 
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Figure 30: Predicted noise level and inputs for a 450 m receiver on 2/27 /91 at 1800. a) 
Noise spectral level 101ogP(/) in dB re yuPa 2 /Hz. b) Wave height power 
spectrum H( a /2 k ) in m 2 /Hz. c) Wave height acoustic input 201ogH(//2) 
in dB re m 2 /Hz. d) Spreading integral level 10 log in dB re Hz” 1 . 
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Figure 31: Meteorological data as measured at Discus E on 1/27/91. a) Wind vector. 

The direction of the arrow indicates the direction of the wind, with a 
downward arrow indicating a wind from the west. The magnitude of the 
arrow indicates the wind speed, with the wind speed scale given on the 
vertical axis, b) Wind speed and wave height. 
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Figure 32: Meteorological data as measured at Discus E on 3/04/91. a) Wind vector. 

The direction of the arrow indicates the direction of the wind, with a 
downward arrow indicating a wind from the west. The magnitude of the 
arrow indicates the wind speed, with the wind speed scale given on the 
vertical axis, b) Wind speed and wave height. 
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Figure 33: Meteorological data as measured at Discus E on 2/22/91-2/23/91. a) 
Wind vector. The direction of the arrow indicates the direction of the 
wind, with a downward arrow indicating a wind from the west. The 
magnitude of the arrow indicates the wind speed, with the wind speed 
scale given on the vertical axis, b) Wind speed and wave height. 
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Figure 34: Amplitude response of the anti-aliasing and pre-whitening filter used in 
the ECONOMEX instruments with an 8 Hz sampling rate. Response 
amplitude in dB re 1 Volt. 
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Figure 35: An example of a full spectral estimate using data from OBS 58 DPG on 
3/04/91 at about 1700 hours. This spectrum was generated as described 
in the text with the exception that a 2048 point FFT was used versus a 
512 point FFT. The spectral level is in dB re /iPa 2 /Hz. 
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Figure 36: Predicted and observed noise lev- 

els, 1/27/91, 1800, OBS 58 (2500 
m), in dB re fi Pa 2 /Hz. Observed 
spectrum is solid curve. 



450 m, 2/23/91, 1400 




Figure 37: Predicted and observed noise lev- 

els, 2/23/91, 1400, OBS 63 (450 
m), in dB re /iPa 2 /Hz. Observed 
spectrum is solid curve. 
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Figure 38: Predicted and observed noise lev- 

els, 2/23/91, 0700, OBS 58 (2500 
m), in dB re /zPa 2 /Hz. Observed 
spectrum is solid curve. 



2500 m, 1/27/91, 1200 




Figure 39: Predicted and observed noise lev- 

els, 1/27/91, 1200, OBS 58 (2500 
m), in dB re /xPa 2 /Hz. Observed 
spectrum is solid curve. 
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Figure 40: Predicted and observed noise lev- 

els, 3/04/91, 0000, OBS 58 (2500 
m), in dB re /iPa 2 /Hz. Observed 
spectrum is solid curve. 
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Figure 41: Predicted and observed noise lev- 

els, 2/27/91, 1800, OBS 63 (450 
m), in dB re /xPa 2 /Hz. Observed 
spectrum is solid curve. 
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Figure 42: Predicted and observed noise lev- 

els, 2/27/91, 1500, OBS 63 (450 
m), in dB re /xPa 2 /Hz. Observed 
spectrum is solid curve. 
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Figure 43: Predicted and observed noise lev- 

els, 2/23/91, 1000, OBS 63 (450 
m), in dB re /xPa 2 /Hz. Observed 
spectrum is solid curve. 
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Figure 44: Predicted and observed noise lev- 

els, 1/27/91, 0300, OBS 58 (2500 
m), in dB re /iPa 2 /Hz. Observed 
spectrum is solid curve. 
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Figure 45: Predicted and observed noise lev- 

els, 1/27/91, 0600, OBS 58 (2500 
m), in dB re /i Pa 2 /Hz. Observed 
spectrum is solid curve. 
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Figure 46: Predicted and observed noise lev- 

els, 1/27/91, 0900, OBS 58 (2500 
m), in dB re /i Pa 2 /Hz. Observed 
spectrum is solid curve. 



Figure 47: Predicted and observed noise lev- 

els, 3/04/91, 0400, OBS 58 (2500 
m), in dB re /iPa 2 /Hz. Observed 
spectrum is solid curve. 
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Figure 48: Predicted and observed noise lev- 

els, 3/04/91, 0600, OBS 58 (2500 
m), in dB re /iPa 2 /Hz. Observed 
spectrum is solid curve. 
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Figure 49: Predicted and observed noise lev- 

els, 3/04/91, 0800, OBS 58 (2500 
m), in dB re /zPa 2 /Hz. Observed 
spectrum is solid curve. 
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Figure 50: Predicted and observed noise lev- 

els, 3/04/91, 1000, OBS 58 (2500 
m), in dB re /*Pa 2 /Hz. Observed 
spectrum is solid curve. 
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Figure 51: Predicted and observed noise lev- 

els, 2/23/91, 0200, OBS 58 (2500 
m), in dB re /^Pa 2 /Hz. Observed 
spectrum is solid curve. 
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Figure 52: Predicted and observed noise lev- 

els, 2/23/91, 0300, OBS 58 (2500 
m), in dB re /iPa 2 /Hz. Observed 
spectrum is solid curve. 
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Figure 53: Predicted and observed noise lev- 

els, 2/23/91, 0400, OBS 58 (2500 
m), in dB re /iPa 2 /Hz. Observed 
spectrum is solid curve. 
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Figure 54: Predicted and observed noise lev- 

els, 2/23/91, 0500, OBS 58 (2500 
m), in dB re /iPa 2 /Hz. Observed 
spectrum is solid curve. 
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Figure 55: Predicted and observed noise lev- 

els, 2/23/91, 0600, OBS 58 (2500 
m), in dB re /iPa 2 /Hz. Observed 
spectrum is solid curve. 



105 



Noise Power (dB) 



2500 m, 1/27/91, 1500 




Frequency (Hz) 



Figure 56: Predicted and observed noise lev- 

els, 1/27/91, 1500, OBS 58 (2500 
m), in dB re /iPa 2 /Hz. Observed 
spectrum is solid curve. 



2500 m, 1/27/91,2100 




Figure 57: Predicted and observed noise lev- 

els, 1/27/91, 2100, OBS 58 (2500 
m), in dB re /iPa 2 /Hz. Observed 
spectrum is solid curve. 
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